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ABSTRACT 

Context. Gravitational lensing is one of the leading tools in understanding the dark side of the Universe. The need 
for accurate, efficient and effective methods which are able to extract this information along with other cosmological 
parameters from cosmic shear data is ever growing. COSEBIs, Complete Orthogonal Sets of E-/B-Integrals, is a recently 
developed statistical measure that encompasses the complete E-/B-mode separable information contained in the shear 
correlation functions measured on a finite angular range. 

Aims. The aim of the present work is to test the properties of this newly developed statistics for a higher-dimensional 
parameter space and to generalize and test it for shear tomography. 

Methods. We use Fisher analysis to study the effectiveness of COSEBIs. We show our results in terms of figure-of-merit 
quantities, based on Fisher matrices. 

Results. We find that a relatively small number of COSEBIs modes is always enough to saturate to the maximum 
information level. This number is always smaller for Jogarithmic COSEBIs' than for Jinear COSEBIs', and also depends 
on the number of redshift bins, the number and choice of cosmological parameters, as well as the survey characteristics. 
Conclusions. COSEBIs provide a very compact way of analyzing cosmic shear data, i.e., all the E-/B-mode separable 
second-order statistical information in the data is reduced to a small number of COSEBIs modes. Furthermore, with 
this method the arbitrariness in data binning is no longer an issue since the COSEBIs modes are discrete. Finally, the 
small number of modes also implies that covariances, and their inverse, are much more conveniently obtainable, e.g., 
from numerical simulations, than for the shear correlation functions themselves. 

Key words. Gravitational lensing- cosmic shear: COSEBIs- methods: statistics 



1. Introduction 

As light travels through the Universe, the gravitational po- 
tential inhomogeneities distort its path; these distortions 
result in sheared galaxy images and carry invaluable in- 
formation about the matter distribution between the ob- 
server and the source. Cosmic shear analysis is the study 
of the effects of large-scale stru ctures on light bundles (see 
iBartelmann fc Schneider] 1200 it ). Consequently, it is one of 
the most promising probes for understanding the Universe, 
especially dark energy. The upcoming cosmic shear surveys 
(e.g. Pan-STARRl© KID93, DE9l LSST0, and Euclicfl) 
will have better statistical precision compared to current 
surveys, which means lower noise levels, larger fields of 
view, deeper images, and more accurate redshift estima- 
tions. Trustworthy and accurate methods are able to ex- 
tract all the potential information in these future observa- 
tions and make the effort put into launching them worth- 
while. 

The most direct second-order statistical measurement 
from any weak lensing survey are the shear two-point cor- 



relation functions £± (■&) , which in reality can be determined 
only on a finite interval # m in < i? < i? max . These, however, 
cannot be used for a comparison with theoretical models, 
since the shear field is in general composed of two modes: 
B-modes cannot be due to leading-order lensing effects, 
although they provide a measure of other effects such as 
sha pe measurement errors and intri nsic alignment effects 



Joachimi fc Schneiderll2010t also ISchneider et al.lfl998l 



Schneider et al.ll2002bl for other effects). On the other 
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(see 
and 

hand, E-modes are the only relevant modes when it comes 
to comparing the cosmic shear data with models. 

Almost all of the recent analysis of cosmic shear 
data employ metho ds of E-/B-mode separation (e.g. 
iBeniamin et alJl2007l and lFu et al.ll2008f ). These studies are 
done in either Fourier or real space. For Fourier space analy- 
sis one has to find an estimate of the power spectrum, which 
is sensitive to gaps and holes in the survey and in general 
the survey geometry, which complicates such analysis. On 
the other hand the studies in real space do not share the 
same complications, since estimators of the shear correla- 
tion functions are unaffected by such gap s. Most of these 
studi es use the aperture mass dispersion (|Schneider et al.l 
Il998l h which applies compe nsated circular fil t ers to the 
shear field. As was show n in ICrittenden et al.l (|2002f ) and 
ISchneider et alj (|2002alh the aperture statistics, in prin- 
ciple, cleanly separates the shear two-point correlations 
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(2PCFs) into E-/B-mode contributions. Furthermore, in 
the two papers just mentioned, a decomposition of the shear 
2PCFs into E- and B-mode correlation function £e,b($) has 
been derived, which also has been employ ed in cosmic shear 
analyses of survey data (|Lin et al.ll2011l i. 

However, both the aperture statistics and the E-/B- 
mode correlation functions are unobservable in practice. 
The aperture mass dispersion requires shape measurements 
of galaxy pairs down to arbitrarily small angular scales. 
Since this is not feasible in real data, usually ray-tracing 
simulations fill i n the gap, resulting in biases and E-/B- 
mode mixing (see lKilbinger et aDl2006T ). On the other hand, 
the determination of £e,b($) requires the knowledge of 
out to infinite Hence, in both cases, determin- 
ing E-/B-mode separated statistics requires some sort of 
data invention. 

To overcome these problems, ISchneider fc Kilbingerl 
(|2007l ) derived general conditions and relations for E- 
/5-statistics based upon two-point statistical quantities, 
namely 2PCFs and convergence power spectra. They de- 
fined the quantities 

E=\ /~dtftf[T4tf)£ + 0?)+T_0?)£_(#)], (1) 



B 



10 



dtftf [T + (<?)£+ (tf)-T_ (<?)£_(<?)]; (2) 



provided that the filter functions satisfy 



d0 0T + (0)Jo(i0) 



d0T?r_(tf)J 4 (&?) , (3) 



E depends only on the E-mode shear, and B depends only 
on the B-mode shear (with the aperture dispersion being 
one particular example). Moreover, they have shown that 
in order to obtain these statistics from the shear 2PCFs on 
a finite angular interval, < i? m i n < i? < # max < oo, the 
filter function T + should have finite support on the same 
angular interval and satisfy 



dtfi?T + (tf) = 



dtftf 3 T + 0?). (4) 



Whereas all solutions to the above relations provide 
statistics which cleanly separate E-/B-modes on a finite 
interval, different solutions may vary in their information 
contents. For example , the ring statistics introduced in 
ISchneider fc Kilbingerl (|2007f ) has a lower signal-to-noise 
for a fixed angular range than the aperture dispersion, 
which, however, is compensated by its more diagonal noise- 
covariance matrix resulting in comparable Fisher ma trices 
with aperture mass dispersion (|Fu fc Kilbingerll2010t) . 

Rec ently, a complete sol ution of this issue was ob- 
tained (jSchneider et al.l 120101 hereafter SEK) by defining 
Complete Orthogonal Sets of E-/B-Integrals (COSEBIs). 
COSEBIs capture the full information of the shear 2PCFs 
on a finite interval which is E-/B-mode separable. In fact, 
SEK have shown that a small number of COSEBIs contain 
all the information about the cosmological dependence in 
their two-parameter model. Furthermore, they showed that 
COSEBIs in fact put tighter constraints on t hese parame - 
ters compared to the aperture mass dispersion. lEiflerl (|201ll ) 
obtained a similar conclusion for a five-parameter cosmolog- 
ical model. Therefore, the set of COSEBIs not only capture 



the full information, but also provide a highly efficient and 
simple method for data compression. 

In this paper we further generalize the analysis in SEK 
to seven cosmological parameters, erg, Q m , f^A, wq, n s , h, 
and 17b, and investigate the effect of tomography on the 
results. Tomography, the joint analysis of shear auto- and 
cross-2PCFs of galaxy populations with different redshift 
distributions, is a powerful tool for cosm ological analysis 
(jAlbrecht et al.ll2006l : IPeacock et alj|2006l). in particular in 
multi -dimensional parameter space (see ISchrabback et ahl 
l20l"0l for a recent paper on constraints on dark energy from 
cosmic shear analysis with tomography). We use Fisher 
analysis throughout our paper to represent the constrain- 
ing power of COSEBIs, and compare the results from a 
medium-sized with that of a large cosmic shear survey. 

In Sect. [5] we summarize the method used in SEK and 
write the corresponding relations for shear tomography. In 
Sect. [3] we briefly explain our choice of cosmology, and in 
Sect.[l]the covariance of COSEBIs is shown. We present our 
figure-of-merit based on Fisher analysis and show the re- 
sults for the seven cosmological parameters and up to eight 
redshift bins in Sect. [5] Finally we conclude by summarizing 
the most important results of the previous sections and em- 
phasizing the advantages of COSEBIs over other methods 
of cosmic shear analysis. We have also derived an analytic 
solution to the linear COSEBIs weight functions presented 
in App.S 

2. COSEBIs 

There is an infinite number of filter functions T+($) satisfy- 
ing Eq. (U) . Such filters can be expanded in sets of orthog- 
onal functions, labeled T+ n (#); the corresponding T_ n (i?) 
are obtain ed from solving Eq. which can be inverted 
explicitly (jSchneider et al. l2002al) . Accordingly, the corre- 
sponding E/B-statistics are denoted by E n and B n , re- 
spectively. Here we will also consider the case that dif- 
ferent galaxy populations can be distinguished (mainly by 
their redshifts); therefore, one can measure auto- and cross- 
correlations functions of the shear, We denote the 
corresponding COSEBIs by E% and B% . They are related 
to the auto- and cross-power spectra of the convergence, by 

dee L 

L E 



ES 



B 13 = 



2n 

dee 



- B (W), 



-P • 



(5) 
(6) 



where P^ B are the E-/B-cross convergence power spectra 

of galaxy populations i and j ( see ISchneider et al"1l2002al ). 
and are related to the 2PCFs by 

de 



em 
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dee 

~2T 



J («)[P E J (£) 



J 4 (W)[i^(*) 



POT 



POT 



(7) 



(8) 



Inserting the above relations into Eq. ([T]), one can find re- 
lations connecting W n to T± n 

W n (e) = / di?7?T +n (0)Jo(&?) (9) 



di?i?r_ n (0)j4(&?) 



(io) 
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Any type of cosmic shear analysis needs some sort of 
error assessment. In particular Fisher analysis, used in the 
present work, depends on the noise-covariance of the statis- 
tics employed. The noise-covariance of COSEBIs for sev- 
eral galaxy populatio ns assuming Gaussian shear fields (see 
Uoachimi etaLll2008[ ) is 



c X(ij,kl) 



{XliX 
1 



kl\ 



2ttA 



-(xn){x%) 

d££W m (£)W n (£) 



where 



P${£) 



Sik 

2n,- 



(11) 



(12) 



and X stands for either E or B. The survey parameters are 
also included in Eq. (fTT|) with the survey area, A, the galaxy 
intrinsic r.m.s ellipticity, <r £ , and the mean number density 
of galaxies in each r edshift bin, n,. 

In a recent paper. lSato et aL (12011) have shown that the 
Gaussian covariance model in Joachimi et al.l (|2008f ) overes- 
timates the true Gaussian covariance for surveys with small 
area (A < lOOOdcg 2 ), and they have developed a fitting for- 
mula to correct for this discrepancy ; in spite of their find - 
ings we will stick to the estimation of lJoachimi et al.l (|2008fh 
since the fitting formula in the latter paper depends on 
source redshift and is developed for a single source galaxy 
redshift, making it non-applicable for this work. 

Alternatively, one can write the covariance (Eq. fTT]) in 
terms of T± n and the two-point correlation functions' co- 
variance (see SEK). However, in this approach double inte- 
grals over the covariance of 2PCFs slow down the calcula- 
tions. 



2.1. The COSEBIs filter and weight functions 

SEK constructed two complete orthogonal sets of functions, 
linear and logarithmic COSEBIs (hereafter Lin- and Log- 
COSEBIs respectively), by considering Eq. (U), and impos- 
ing orthogonality conditions on the T +n filters. Once the 
T +n filters are known, the T_„ filters can be calculated via 
Eq. ©. The Lin-COSEBIs filters are polynomials in 0, the 
angular separation of galaxies, while the Log-COSEBIs fil- 
ters are polynomials in ln(i9). 

The output of theoretical cosmological models which is 
of relevence here is the power spectrum. Hence, the quick- 
est way to treat COSEBIs in theory is to work in ^-space 
and to use Eq. (jTTJ) for the covariance, without taking the 
detour of calculating the shear 2PCFs and their covariance. 
As a result we need to calculate the W n (£) functions which 
are the Hankel transform (Eq.|H]) of their real-space counter- 
parts, T± n . For convenience, we choose to evaluate W n (£) 
from their integral relation with Jo and T +n . Since both 
Jo and T +n are oscillating functions, evaluating these inte- 
grals is rather challenging, in particular for large £. A piece- 
wise integration, from one extremum to the next, is used in 
the present work to evaluate W n (£). App.lAl contains more 
details about the numerical integrations and also a (semi- 
analytic formula for the linear W n functions. 

As is explained in SEK, the Log-COSEBIs are more ef- 
ficient for a cosmic shear analysis. The reason is that unlike 



Table 1. The fiducial cosmological parameters consistent with 
WMAP 7-years results. The normalization of the power spec- 
trum, ag, is the standard deviation of perturbations in a sphere 
of radius 8 fe~ 1 Mpc today. fi m , Qa, and fib are the matter, 
the dark energy and the baryonic matter density parameters, 
respectively, wo is the dark energy equation of state parameter, 
which is equal to the ratio of dark energy pressure to its den- 
sity. The spectral index, n s , is the slope of the primordial power 
spectrum. The dimensionless Hubble constant, in Ho — lOOh 
km s _1 Mpc -1 , characterizes the rate of expansion today. 
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the linear filter functions which oscillate fairly uniformly in 
linear scale, the logarithmic T+° s have their roots fairly uni- 
formly distributed in log($), i.e., they are more sensitive to 
variations of £± on smaller scales. Combining this property 
with the fact that most of the cosmic shear information 
is contained in these smaller scales shows that it is more 
reasonable to employ Log-COSEBIs. In the next section 
we will show the difference of the Log- and Lin-COSEBIs 
using our figure-of-merit. 

In Fig. Q] and Fig. [2] the behavior of linear and logarith- 
mic COSEBIs weight functions, W^ in (£) and W^° s (£), for 
three angular ranges can be seen. The W^° e (£) and W^ ln (£) 
have different yet similar oscillatory properties. They both 
die out rapidly with increasing £ but the lower frequency 
oscillations of W^° s (£) are more prominent. They show ap- 
proximately the same inverse relation to $ max and $ m i n for 
their lower and upper limits. 

3. Cosmological Model 

The cosmological model assumed in the present work is 
a wCDM model (jPeebles fc Ratral 120031 and references 
therein), i.e., a cold dark matter model including a dynam- 
ical dark energy with an equation-of-state parameter wq. 
The fiducial value of the parameters involved are listed in 
Tab.m 

The starting point in the analysis is to derive the matter 
power spectrum. Fo r the li near power spectrum we used the 
iBond fc Efstathioul (1 19841) trans fer function, and the halo 
fit formula of ISmith et all (|2003h for a fit of the non-linear 
regime. 

To calculate the convergence power spectrum we need 
the redshift distribution of galaxies. The overall redshift 
probability distribution is parametrized by 



p(z) 







2 T[(l + a)//3] 



cxp 



(13) 



which represents the galaxy distribut ion fairly well (it is a 
generalization of lBrainerd et al.1ll996t ). The parameters, a, 
/3, and zq depend on the survey. We consider a medium and 
a large survey (hereafter MS and LS respectively). The MS 
has the same area as the CFHTLS (|Fu et al.ll2008l ). a cur- 
rent survey, while the LS covers the whole extragalactic sky 
and represents future surveys. The parameters of our two 
model surveys are given in Tab.[5J and the corresponding 
redshift distributions are plotted in Fig.[3j 

Constructing the Fisher matrix requires the derivatives 
of the E-mode COSEBIs and of their covariances with re- 
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Fig. 1. The weight functions Wn ln (£) are the Hankel transforms of T± n (i}) as in Eq. (J9j . In the blow-ups, the two modes of 
oscillation for each W^ ln can be seen, the lower frequency mode and the higher frequency mode which are inversely proportional 
to i? m in and $max, respectively. The overall amplitude of the oscillations strongly depends on n and i9 ma x- 
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Fig. 2. The weight functions Wn° s (£) are the Hankel transformation of T^ og (tf) as in Eq. ©. Similar to the W^ in , the position 
of the first peak depends mainly on i9 max and is rather insensitive to $ m i n - The difference between the two sets of linear and 
logarithmic function can be seen most prominently in the blow-ups; the lower frequency oscillations are more pronounced in this 
case. 
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Table 2. The redshift distribution parameters and the survey parameters for our medium and large surveys, a, (3, and zo determine 
the total redshift distribution of sources, while z m i n and 2 max indicate the minimum and the maximum redshifts of the sources 
considered. A is the survey area in units of deg 2 , a € is the galaxy intrinsic ellipticity dispersion, and n is the mean number density 
of sources per square arcminute in the field. 
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Fig. 3. The overall source redshift probability distribution of 
source galaxies assumed for the two surveys. LS has a deeper 
source distribution compared to MS. 



spect to the parameters. For example, to take the derivative 
with respect to fi m , its relation to the shape par ameter, T 
should be notified. In the present work we use the lSugivamal 
(|1995f) relation, 



T = fi m /iexp[-rj b (l + V2h/tt m )\ . 



(14) 



In their derivatives with respect to f2 m , SEK assumed a con- 
stant r, equivalent to allowing h or J7b to vary accordingly 
(the only dependence of the convergence power spectrum 
on h or fib comes through T). In the present work h and 
fib are independent parameters and T depends explicitly 
on f2 m . The difference between the two approaches is not 
negligible, as shown in Fig. 2] which displays the derivative 
of the power spectrum with respect to f2 m in both cases. 
This difference is due to the non-linear relation between 
h and T. To justify our choice of parametrization, we just 
mention that the constraints from cosmological probes on h 
is tighter compared to T, and that makes it a more natural 
choice especially when priors are used. 



4. COSEBIs Covariance 

Fig. [5] shows the noise-covariance of linear and logarithmic 
E-mode COSEBIs for the model parameters of the MS (the 
covariance has a similar behavior in the case of the LS but 
with a different amplitude). This covariance is calculated 
from Eq. (jTTJ) assuming a single source redshift distribution 
(Ea. rnt)) . Moreover, by defining the correlation coefficients 



Fig. 4. The absolute value of the derivative of the convergence 
power spectrum with respect to f2 m . Both of the curves rely on 
a five point stencil method where 4 nearby points have to be 
evaluated. The solid curve is drawn assuming all parameters are 
fixed except Sl m and T, in contrast to the dotted curve where 
instead of V, h or Qh are variable. 



of COSEBIs, 



tain 



AIN 



AIM NN 



(15) 



the behavior of the off-diagonal terms becomes clearer. 
(The capital subscripts N and M can be different from 
the COSEBIs subscripts, if several source populations are 
considered; see below for more details.) Fig. [6] compares the 
correlation coefficients for three different choices of the an- 
gular range, [l',400'], [20', 400'], and [1',20'j, at a fixed 
M = 9. 

Cosmic shear analysis, as we will see in Sect.[5J provides 
more information when redshift information is available. In 
practice the redshifts of galaxies are estimated using sev - 
eral photometric filters (see e.g. iHildebrandt et aT]|2010t ), 
from which an overall distribution for the source galaxies 
is obtained. The distribution is then divided into a number 
of photometric redshift bins. The photometric redshifts are 
not exact, so the true redshift distributions will overlap. 
Therefore, instead of redshift bins, in general one has to 
consider redshift distributions. However, for simplicity, in 
the present work we have assumed redshift bins with sharp 
cuts and no overlap. In addition, the bins are selected such 
that the number of galaxies in each bin is the same. 

In general, a tomographic covariance for r redshift bins 
consists of [r(r + l)/2] 2 building blocks, each of which is a 
covariance matrix of (E% , E^ ) where i, j, k, I are fixed and 



1,2, 



This means in total the covariance 
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Fig. 5. A 3D representation of the non-tomographic covariance 
of 15 E-mode COSEBIs for an angular range of [1', 400'], for MS 
parameters. The x- and y- axes correspond to the elements of the 
covariance matrix, and the value of the vertical axis shows the 
value of the covariance of the corresponding element. A contour 
representation of the covariance is shown for each plot at its 
base. 



matrix has [r(r + 1) 7i m ax] 2 /4 elements, where n max is the 
maximum number of COSEBIs modes considered. 

Nevertheless, a covariance matrix is by definition sym- 
metric and a tomographic covariance is made up of smaller 
covariances, i.e., only x(x + l)/2 x rt max (n max + l)/2 ele- 
ments, with x = r(r + l)/2, have to be calculated, the rest 
are equal to these (see Fig.0). 

The covariance of the E% depends on six indices; in 
order to apply normal matrix operations, the three indices 
of are combined into one 'superindex' N, given by 



N 



(i — 1) x r- 



(i-m -2) 



+u-i) 



(16) 



where r is the total number of redshift bins and n max is the 
total number of COSEBIs modes. 

Using the new labeling, the correlation coefficients of 
E} 1 and i?| 3 (corresponding to N = 7 and N = 82, re- 
spectively) with the other E^ is shown in Fig. [5J where 15 
COSEBIs modes and 4 redshift bins are considered. Each of 
the peaks in the figure correspond to the correlation coef- 
ficient of E] 1 and E!j? . The highest peak with r = 1 occurs 
for M = N, while the rest of the peaks are correlations 
between different redshift bins. The Log-COSEBIs show 



Fig. 6. The correlation coefficients of non-tomographic 
COSEBIs for different angular ranges [t) m i n ,iJmax] at m = 9, for 
the MS parameters. Here M, the capital subscripts, are equal 
to the COSEBIs mode, m. 



larger noise-correlations between different modes compared 
to Lin-COSEBIs, which may persuade one to choose the 
Lin-COSEBIs for cosmic shear analysis. However, the Log- 
COSEBIs compensate this apparent disadvantage by re- 
quiring fewer modes to saturate the Fisher information level 
for relevant cosmological parameters compared to the linear 
ones, i.e., the number of covariance elements that have to 
be calculated for Lin-COSEBIs is higher and hence deter- 
mining their covariance matrix is more time consuming, es- 
pecially when redshift binning is considered. Consequently, 
in Sect. [5] we mainly employ Log-COSEBIs to analyze to- 
mographic Fisher information. 

5. Results of Fisher analysis 

5.1. Figure-of- merit 

In this section we carry out a figure-of-merit analysis to 
demonstrate the capability of COSEBIs to constrain cos- 
mological parameters from cosmic shear data. Our figure- 
of-merit, /, based on the Fisher matrix, quantifies the 
credibility of the estimated parameters. In general, for 
any unbiased estimator, the Fisher matrix gives the lower 
limit of the errors on par ameter estimations (se e e.g. 
iKennev fc Keeping|[l95ll and Kendall fc StTIaItl[T960l for de- 
tails). 
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Fig. 7. A representation of a tomographic covariance. In this di- 
agram 3 redshift-bins (1,2,3) and 5 COSEBIs modes are assumed 
to be present. The blow-up shows one of the covariance building 
blocks; the numbers 1-5 show the COSEBIs mode considered, 
e.g. 15 means the covariance of E\ and E$. The numbers on the 
sides of the matrix show which combination of redshift bins is 
considered, e.g., 12 means the covariance of redshift-bins 1 and 
2 is relevant. Due to symmetry, only a part of the covariance 
elements have to be calculated, here shown in pink. 
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The Fisher matrix is related to the COSEBIs by 

Fij = \^[C- 1 Ci C- 1 C d + C- 1 M, l3 ] , (17) 
where C is the COSEBIs 

E. } Efi, E is the vector of the E-mode COSEBIs, and 
the commas followed by subscripts indicate partial deriva- 
tives with respect to the cosmological parameters (see 
iTeemark et al.lll997l for example) . We define our figure-of- 
merit, /, in a very similar manner to SEK 



/ 



1 



Vdct F 



l/r. 



where n p is the number of free parameters considered. In 
the following analysis, we will assume for simplicity that 
the first term in Eq. (]17p is much smaller than the second 
and can thus be neglected. Note that this approximation 
becomes more realistic in the case of a large survey area, 
since the first term on the r.h.s. of Eq. (fTT)) does not depend 
on the survey area, while the second term is proportional 
to it (recall that C oc \/A or in other words C _1 oc A). We 
checked that our medium survey is already big enough for 
this approximation to hold (see Fig. [5]). 

With the definition (fT5]) we compress the Fisher matrix 
into a one-dimensional quantity, which provides a measure 
of the geometric mean of the standard deviations of the 
parameters; e.g. in the case of one free parameter <fi, f{4>) 
is equal to the standard deviation a((f>) of that parameterjj 



6 Another quantity, q, was also defined in SEK to measure the 
area of the likelihood regions. It is calculated from the second- 



Fig. 8. The correlation coefficients of COSEBIs for an angular 
range of [l',400'] and 4 redshift bins. In total, 15 COSEBIs 
modes are considered for each graph. The rjuiv is shown for 
M — 7 corresponding to Ej 1 , and for M = 82 corresponding to 
Ef. 
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Fig. 9. A comparison between a simplified and complete Fisher 
analysis, using Log-COSEBIs. The asterisks show the case where 
the derivatives of the covariance is taken into account (the first 
part of Eq . I17fl while the squares show the simplified case where 
we assume these derivatives are zero, in calculating /. Here as 
is the only free parameter, whereas the rest of the parameters is 
fixed to their fiducial values. 
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Fig. 10. A comparison between the Log- and Lin- 
COSEBIs Fisher analysis results for two sets of 
assumptions, where as and Jim are the free pa- 
rameters and the rest is fixed to their fiducial 
values. In one case the shape parameter P is held 
fixed, while in the other it is left as a variable de- 
pending on Q m and the fiducial values of h and Ob 
(according to Eq. ll4[) . The same analysis is also 
carried out for the Full-COSEBIs and the shear 
2PCFs. 
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Fig. 11. These plots show one of our consistency checks. It is a comparison between the Lin- and Log-COSEBIs results for LS 
with a single (left panel) and two galaxy redshift distributions (right panel), including all of the 7 parameters. Apart from very 
small numerical inaccuracies, both sets of COSEBIs saturate to the same value, as expected. There are two solid lines in each plot. 
The line with the higher value shows the value of Log-COSEBIs at n max = 20, and the other line shows the value of / as obtained 
from the shear 2PCFs. The slightly smaller value of / in the latter case (this difference is not visible in the plot) is related to the 
fact that the analysis from the shear 2PCFs implicitly assume the absence of B-modes, and thus contains information from very 
large-scale modes which, however, cannot be uniquely assigned to either E- or B-modes. The comparison of the two plots shows 
that dividing the galaxies into two redshift bins not only increases the information content of the Fisher analysis but also decreases 
the number of COSEBIs modes needed. Note that the x-axis of the single redshift distribution plot starts from 7, the other one 
from 3. 



5.2. Assumptions and parameter settings 

For our cosmic shear analysis we considered a medium (MS) 
and a large survey (LS) as explained in Sect. [3] We have 
also studied the effect of a Gaussian prior, in the form of 
a Fisher matrix. This prior is the inverse of the WMAP7 
parameter covariance matrix from the final iteration (5000 
sam ple points) of a Popu lation Monte Carlo (PMC) run 
(see lKilbinger et aL1l2010l ). called the CMB prior from here 



order moments of the posterior likeliho od, q and / are equal if 
the posterior is a multivariate Gaussian. lEiflerl (|20 1 if ) has shown 
that the difference between / and q is small, especially for a large 
survey area. 



on. We implement the prior by adding the Fisher matrices 
of our COSEBIs analysis and the CMB prior. The value of 
the CMB prior is shown in terms of f((f>) in the first column 
of Tab. [3] for each of the parameters. 

We consider three different angular ranges, [1', 20'] , 
[l',400 '] , and [20', 400'] . The motivation for this choice 
is as follows: We consider a total interval of [1', 400'] where 
the flat sky approximation is still valid up to the maxi- 
mum separation and galaxy shapes are easily distinguish- 
able for the minimum separation; also, i9 m i n = 1' avoids 
the scales where baryonic effects are expected to have the 
strongest ifluence. We further divide this interval into two 
non-overlapping parts with $ m ax/$min = 20, to compare 
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cosmic shear information on small and large scales. The 
small-scale range, [1', 20'] may apply for a cosmic shear 
survey of individual one square degree fields. The large 
scale interval. [20', 400'] could be used for very conserva- 
tive analyses where non-linear and baryonic effects are to 
be avoided. 

In Sect. 15.31 we show the value of / for two parameters 
while the rest are fixed to their fiducial values for the MS, 
and also for all seven parameters for the LS. In principle we 
could show all of the possible combinations for parameters, 
nevertheless finding the error on each of the parameters 
seems a more relevant task. Therefore, the rest of our anal- 
ysis, carried out in Sect. l5.4l is done for a single parameter, 
4>, where f{<j>) = 0$. 

To find the value of / for a single parameter we use two 
approaches. In one approach we fix the six other parameters 
to their fiducial values in Tab.[JJ while in the other case we 
marginalized over the remaining six parameters. 

For each setup we investigate the amount of information 
with respect to the number of COSEBIs modes considered. 
In addition we analyze the behavior of / with the number 
of redshift bins considered. 

5.3. Properties of COSEBIs 

For a fixed number of modes, the Log-COSEBIs are more 
sensitive than the Lin-COSEBIs to structures of the shear 
2PCFs on small scales. Here we show its effect on the Fisher 
analysis. We investigate the dependence of / on the number 
Umax of COSEBIs modes incorporated in the analysis. 

SEK have shown the difference between the behavior 
of the Lin- and Log-COSEBIs for two parameters, eg and 
Q m , with T fixed (their definition of / and fiducial values 
of parameters are slightly different from ours). Similar to 
their work, we here compare the values of / for the same 
two parameters with Lin- and Log-COSEBIs. In addition, 
we inspect the difference between a fixed shape parameter, 
or its dependence as given in Eq. ([T4"]). 

Fig. [10] is a representation of our inspection for the MS 
in the angular range of [l',400'] . Two general conclusions 
come out of this comparison: (1) / for fixed and depen- 
dent r converges to the same value for the Lin- and Log- 
COSEBIs. (2) The values of / for a fixed or dependent T are 
different, and also the convergence rate is different. E.g., the 
Lin-COSEBIs reach the saturated / value for n max ps 40 
for a variable T, while in the other case, only 25 modes 
are needed. This effect is less dramatic in the case of Log- 
COSEBIs (they need 7 modes for a variable V and 5 modes 
for a constant one), since they generally converge faster. 
Similarly in Fig.[TTJ we visualize our consistency check by 
showing that the values of / for Log- and Lin-COSEBIs 
converge to the same value for seven parameters 

In Figs.[TTj]and[TTl we also show the value of / as derived 
directly from the shear 2PCFs, i.e., without E-/B-mode sep- 

7 In general there are slight differences between the final value 
of / due to numerical inaccuracies, but these differences never 
exceed a few percent and are typically much smaller. An excep- 
tion happens when the saturation is too slow, and the Fisher 
matrix elements are too small, which is the case for MS with 
one redshift bin and 7 parameters, observable especially after 
marginalizing over 6 parameters when the remaining parameter 
or ag. However, for these cases, / is much larger 
than unity, i.e., cases in which no meaningful constraints can be 
obtained anyway. 



aration. As expected, in this case / becomes slightly smaller 
since it is now implicitly assumed that all the signal is due 
to E-modes. However, this is not justified in general; for 
example, very large-scale modes (i.e., small £) enter £+(1?) 
even for small i9, and such modes cannot be uniquely as- 
signed to either E- or B-modes. Thus, the decrease of /, 
and accordingly, the information gain is just an apparent 
one, bought by making a strong assumption. The relative 
difference between the 2PCFs and the converged Lin-/Log- 
COSEBIs values for / is larger for the variable T case, since 
here small-l? modes, which are filtered out in the COSEBIs, 
contain information about the power spectrum shape. 

We also considered as further possibility that the re- 
quirement of finite support for the £_ ($) is dropped, and 
call this 'Full-COSEBIs'. They form a complete set of 
functions on [$ m in,$max], without the constraints given 
in Eq. Though not physically reasonable, the Full- 
COSEBIs are equivalent to measuring £ + only, on the same 
interval. As can be seen from Fig.[JJjl the full COSEBIs yield 
a slightly lower value of / than the true COSEBIs, showing 
that £_ on scales larger than i9 max adds apparent informa- 
tion, which, however, is not observable. We stress here that 
the E-/B-mode correlation f unctions £e/b> in troduced by 
ICrittenden et alj (|2002T ) and ISchneider et all (|2002al) . are 
essentially equivalent to the Full-COSEBIs, since they are 
also based on the assumption that can be measured to 
arbitrarily large separations - which, however, is not pos- 
sible. Therefore, a cosmic shea r analysis based on £e (e.g., 
iFu et all l2008t iLin et aD 1201 lh underestimates the uncer- 
tainties of cosmological parameters. 

Furthermore, we compare the Lin- and Log-COSEBIs 
for LS parameters in Fig. llll for one and two redshift bins. 
The x-axis in the left plot starts from 7 in contrast to the 
right one which starts from 3. The reason is that to con- 
strain rip parameters at least n p equations are needed, i.e., if 
one redshift bin is considered, n p COSEBIs modes should 
be accounted for to produce a covariance matrix with at 
least n p x n p elements. For more than one redshift bin, 
a smaller number of COSEBIs modes are sufficient, subse- 
quently the saturation rate of / is faster, as is visible in the 
right plot in the figure. Recall that 2 redshift bins means 
3 different redshift combinations, i.e., for 7 parameters, the 
smallest integer not less than [7/3] = 3 COSEBIs modes 
are needed. 

5.4. Forecast for parameter constraints 

This section is dedicated to our final results according to 
the assumptions and parameters explained in Sect. 15.21 

Fig.QJ shows the dependence of / for 20 Log-COSEBIs 
modes and for the [l',400'] angular range on the number 
of galaxy distributions (i.e., redshift bins), where all but 
one parameter are marginalized over. Dividing the galaxy 
distribution into more than 4 redshift bins does not change 
the value of / considerably. Nevertheless, a much larger 
number of redshift bins is required to control and correct 
for systematic e ffects, e.g., coming from intrin sic alignments 
(see for example iJoachimi fc Sc hneider 2 0101 and references 
therein). 



8 They are obtained by adding two additional weight functions 
T+ to those used in the COSEBIs; for the linear case, we just 
take all Legendre polynomials (see SEK). 
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Fig. 12. The dependence of the COSEBIs saturated information content in form of / for each of the seven parameters while the 
rest is marginalized over. The above plots are relevant to MS and LS parameters for an angular ranges of [l',400'] . The redshift 
dependence of / is very outstanding here and hence we needed to use logarithmic scales for the y-axes, especially in the cases 
where cosmic shear analysis is done without any priors. 



We also show the dependence of / on n max , for 8 redshift 
bins and marginalized parameters, in Fig. 1131 Comparing 
the cosmic shear analysis with and without CMB prior, we 
see from the figure that the prior in general flattens the 
curves. However, the curves are flatter for MS+CMB than 
LS+CMB as a result of the larger difference between the 
LS and the CMB prior. 

The constraints on each of the cosmological parameters 
behave differently with respect to the number of COSEBIs 
modes or redshift bins considered. For marginalized param- 
eters where the behavior of parameters is entangled, their 
curves show a similar decline. 

By comparing the different angular ranges we conclude 
that a wider angular range needs more modes to extract all 
information. We also note that the behavior of the seven 
parameters are not similar and each of them should be fol- 
lowed separately. 

Based on the results from these two figures, we will re- 
port additional results for n max = 20, where the value of / 
is converged, and for either one or eight redshift bins. These 
results are shown in Tab.[3]in the form of f( 4>) for different 
cases. We have compared these values with iDebono et al.l 
(|2010f) . and found them fully consistent. 

In the following we explain our conclusions from the two 
mentioned figures and Tab. [3] in more detail: 

— MS vs. LS vs. CMB prior: In general, because of its 
much larger survey area and larger galaxy number den- 
sity, LS puts tighter constraints on all of the parameters 
than the MS. Furthermore, since the LS is deeper than 
the MS, it allows more sensitive constraints on param- 
eters which are sensitive to the growth of structure, in 
particular wq. As can be seen from Fig.[T3l the requested 
number of COSEBIs for saturation is slightly higher for 
the LS since this survey contains more information, but 
smaller than 20 in all cases. 

The MS constraints on parameters are weaker than the 
CMB alone for almost all cases. f2 m and cr^, the two pa- 



rameters for which present cosmic shear studies provide 
the most relevant constraints, are the only two for which 
the MS constraints are comparable with CMB. As for 
the rest of the parameters, except for f(wo, [1', 400']) for 
the case of fixed parameters, the parameter uncertainty 
of the CMB prior is about one order of magnitude or 
more smaller than that of the MS. We conclude that the 
MS is not large enough to be competitive with the CMB 
for constraining a seven parameter cosmological model. 
However, the combination of the two slightly tightens 
the constraints. 

In contrast to the MS, the LS yields parameter con- 
straints which can be considerably stronger than the 
CMB alone, in particular when tomography is em- 
ployed. To wit, for the total angular range of [l',400'] 
and 8 z-bins, the LS errors are smaller than those from 
the CMB, except for n s with marginalized parameters. 
For fixed parameters, the / value resulting from the 
combination of the LS and the CMB prior is very close 
to the LS value, whereas it can be much smaller than 
that of the CMB alone. I.e., we conclude that the re- 
sulting constraints from the LS are not dominated by 
the assumed prior. 

In contrast to the MS, the LS is able to put useful con- 
straints on the parameters, even without tomography 
and for marginalized parameters. This is seen in the 
difference between the error values obtained from CMB 
alone and LS+CMB. 

The relative value of errors on parameters is different 
between the two surveys. The reason is that the redshift 
distributions of the two surveys are different (see Tab. [2] 
and Fig.[3|) as seen in the figures. Using redshift informa- 
tion in general is equivalent to using structure evolution 
information. The large-scale evolution is more visible in 
the case of a wider redshift distribution which starts 
from z = 0, where these structures are more evolved. 
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Table 3. Listed below are the values of / for 20 Log-COSEBIs modes, where the saturation level is well reached for 
all the cases. The first column indicates which parameter is free, the rest of the parameters are either fixed (4th- 7th 
columns) or marginalized over (8th-llth columns). The two numbers written in the first-column box of each parameter 
in parenthesis are the corresponding / values from CMB alone, for fixed and marginalized parameters, respectively. The 
second and third columns show the survey and angular range considered. The values of / for one and eight redshift bins 
are presented here for comparison. MS and LS stand for a medium and large survey, respectively. 
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— Fixed vs. Marginalized: The difference between the 
value of / for fixed and marginalized parameters is im- 
mense, especially when prior information is not avail- 
able. However, we can state that, for all cases without 
priors, Q m is the best constrained parameter. The rel- 
ative value of the parameters is also different between 
fixed and marginalized cases. Also for the fixed parame- 
ter case, the convergence rate of LS is slower compared 
to the MS and its relative information content is higher, 
as expected. 



In contrast to the fixed parameters case, tomography 
substantially lowers the errors for marginalized param- 
eters. 
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-1.4 -1.2 -1.0 -0.8 -0.6 

Fig. 14. Contour plots of la constraints for pairs of parameters. 
The plot legends are abbreviated to save space. M stands for 
Marginalized, C for the CMB prior, and #z for the number of 
redshift bins. Top plot: The la contours of too and Qa for the 
LS. The two smaller contours are for the case where the five 
other parameters are held fixed to their fiducial values. Bottom 
plot: The la contours of as and wq for the MS. 



— Angular ranges: From the table, we find the following 
common trends for fixed parameters: 

f(l', 400') < f(l', 20') < /(20', 400') , (19) 

V /(!', 400') J MS Ul',400')As ' 
/ /(20',400') \ / /(20',400') \ 
U(l'.400') J us U(l',400')As ' 1 ' 

The first relation shows that there is more information 
at smaller scales than in the large-scale angular inter- 
val, although there is some independent information at 
larger scales. More interesting are the second and third 
relations which show that there is more information at 
smaller scales regarding LS compared to MS, which is a 
consequence of their different redshift distributions and 
cuts (see Fig. [3] and Tab. [2]). The relations between the 
/(''min,''max) for different angular ranges change when 
parameters are marginalized over. In this case the above 
inequalities are no longer valid for all of the parameters. 

Fig. H3] shows constraints for two pairs of parameters. In 
the bottom plot the direction of the contours for the MS 
are determined by the CMB prior. In this case tomography 
slightly improves the constraints on both parameters. The 
top plot, on the other hand, shows the constraints on dark 
energy parameters. In this case tomographic improvements 
are more visible. Also here the direction of the contours are 
different for the two cases. 



6. Summary and Conclusion 

We have generalized the cosmic shear analysis with 
COSEBIs to include seven cosmological parameters, and 
investigated the effect of tomography on parameter con- 
straints. For our analysis we mainly used the Log-COSEBIs, 
although for consistency checks we have shown that Lin- 
COSEBIs results agree with their logarithmic counterparts. 
In App.|X]we show that there are analytic solutions for lin- 
ear W n functions, which are in principle less computation- 
ally demanding. 

Besides the number of parameters and the use of red- 
shift information, the main difference to SEK is a technical 
one: we calculated the COSEBIs and their covariance di- 
rectly from the power spectrum, without using the 2PCFs 
as intermediate step. This choice is more convenient for 
theoretical considerations and highly speeds up the calcu- 
lations of the COSEBIs covariances, an advantage in par- 
ticular for the case of tomography. We also confirmed that 
our method can reproduce the results of SEK. 

We investigated the effect of a Gaussian prior on two 
mock surveys, using three angular ranges on which the 
shear 2PCFs are asumed to be measured, by Fisher analysis 
methods. We considered the case that all but one parame- 
ter are fixed, as well as that where we marginalize over the 
other six parameters, in order to find the constraints on a 
single parameter. The prior was the Fisher matrix resulting 
from a Population Monte Carlo (PMC) analysis of WMAP7 
results. We considered a medium and a large mock survey 
resembling the state-of-the-art in cosmic shear and that of 
future all-(extragalactic) sky surveys. 

Most importantly, we found that a relatively small num- 
ber of COSEBIs captures essentially all cosmological infor- 
mation from a cosmic shear survey. Whereas this number 
is larger than in SEK, due to the higher-dimensional pa- 
rameter space and the tomographic analysis, COSEBIs not 
only act as a clean E-/B-mode separating shear statistics, 
but also as a highly efficient data compression method. We 
stress that this feature is extremely useful also for evaluat- 
ing covariances from numerical simulations. 

The required number of COSEBIs to saturate the 
cosmological information is considerably smaller for Log- 
COSEBIs than for Lin-COSEBIs, which implies a clear 
preference for the former. It also increases with the number 
of free parameters, and depends on the parameters consid- 
ered, the survey, and the angular range on which the 2PCFs 
are measured. In all cases we considered, fewer than 20 Log- 
COSEBIs modes were sufficient to reach information satu- 
ration. Aside from the tighter constraints of the large sur- 
vey (LS) on all of the parameters compared to the medium- 
sized survey (MS), the order of the parameters with respect 
to their Fisher information is different between the two 
surveys. Moreover, in general LS requires more COSEBIs 
modes due to its higher information level. The comparison 
of the three angular ranges shows that most of the informa- 
tion in cosmic shear is contained at smaller scales which is 
why the Log-COSEBIs with their finer oscillations towards 
smaller scales are more sensitive and reach the saturated 
level of information with fewer modes. However, there is in- 
teresting independent information at larger scales, resulting 
in tighter constraints when using both angular ranges. 

We have investigated the dependence of our figure of 
merit. /, on the number of redshift bins considered. In 
agreement with earlier work, we found that tomography 
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greatly tightens parameter constraints, but the cosmologi- 
cal information saturates at around three or four redshift 
bins. This, however, does not imply that coarse redshift 
information is sufficient for future lensing surveys, since 
good redshift information is required to eliminate system- 
atica from the data, such a s intrinsic alignment effects 
fe.g..lKing fc Schneiderll2003l : I.Toachimi fc Schneiderll2008t 
Uoachimi fc Bridldl201oF ~ 

For future work it will be interesting to investi- 
gate t he effects of nulling with C OSEBIs. Nulling tech- 
niques (|Joachimi fc Sch neider 2008) eliminate the intrinsic- 
intrinsic and intrinsic-shear correlations from observed el- 
lipticity correlations. The intrinsic-intrinsic correlation can 
be handled by accurate redshift information to eliminate 
pairs with physical connections. Consequently, one can in- 
vestigate how the cosmic shear information evolves by do- 
ing so, and how many redshift bins are needed in this 
case. Furthermore, our assumption of a Gaussian covari- 
ance becomes unrealistic at small angular scales; hence, it 
will be interesting to carry out a similar analysis based 
on more realistic covariances, either obtained from ray- 
tracing through cosmological density fields or using (semi- 
) analytic models, such as based on log-normal fields (see 
iHilbert et al1l201lh . 
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Appendix A: Analytic solutions to linear COSEBIs 
weight functions 

Calculating the W n functions and evaluating integrals in- 
volving them requires careful methods, as a result of their 
very oscillating nature. In this section we first show our 
semi-analytical solutions to W^ ln (£) and at the end discuss 
our method of integration in more detail. 

The filters, T|jf, were calculated in SEK. By a simple 



variable change of y — M, Eq. ([9J becomes 



dyy T +n{y/t)2o(y) 



One can write T +n in the form 

n+l 



T+„(i9) = ^a m ^, 



(A.l) 



(A.2) 



n+l 



2=0 



)-Si(M min )] 



(A.5) 



The S n functions can be obtained using standard Bessel 
functions relations, which in our case specialize to 



dy 



[y Ji(y)] = y J o(y) 



J- 1 (y)-Ji(»)=2J{,(y) 

J-i(y) = -Mv) 
-Mv) = My) 



(A.6) 

(A.7) 
(A.8) 



where the last equation results from the two equations be- 
fore it. 

Using Eqs. (|A.6[) and (|A.8|) . the following calculations 
for S n are carried out, 



S n (b)= / dyy n+1 3 (y) 
Jo 

= y n+1 Mv) b ~n I &xy n 3i{y) 



= y n+l My) +ny n My) ~n 2 dyy n - L J (y), 
o Jo 

(A.9) 

where we have carried out integration by parts twice. Note 
that the last term in the above equations is equal to S n ^2- 
Consequently, the recursive formula for S„ is 



S n (b) = b n+1 J 1 (b) + n b n J {b) - n 2 5„_ 2 (6) . (A. 10) 
The first two of these functions are 



S (b)= / dyy.h(y)=b.h(b) 



(ATI) 



and then rewrite W„ as: 



SLi(6)= / d2/J (?/) = 6iF 2 [lAl,3/2,-6 2 /4], (A.12) 
Jo 

where iF? is the Hypergeometric functions (see 
lArfken fc Weberl Il995f ). Although there is an analytic 
formula for S—i, it is more convenient to solve it numeri- 
cally using a similar step-wise integration method as was 
explained in Sect. 12.11 . The only difference here is that the 
steps are taken between zeros of the Bessel function. A 
plot of 5_i can be seen in Fig. lAII 

The recursive formula in Eq. ([A.lOp can be rewritten as 
a closed form formula of the form: 



■ni 



i=0 



dy y l+i Mx) 



We define the functions S n as 



S n (b) := / dxy n+1 3 a (y) 



Inserting the above equation into Eq. (|A.3jl gives 



(A.4) 



N k-1 



(A.3) S n (b) = ,h(b)[b n+1 +J2(-^) k Y[(n-^) 2 b n+1 ' 2k } 

k=l i=0 
N k-1 

+ Mb)[nb n + E(-!) fe U( n ~ 2i ) 2 ( n - 2fc)6"- 2fc ] 

fe=l i=0 

N 

+ (-1) N+1 U(n - 2zfS n _ 2(N+l) {b) , (A.I3) 



i=0 
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1000 



Fig. A.l. The weight functions W„(£) depend on Sn 

integrals, 

which are recursively related to each other. Here the shape of 
the starting function S-i is shown. As is expected it reaches 
zero as the argument gets smaller and has a highly oscillating 
behavior similar to its integrand, Jo. 



where N is 



N 



n- 2 
2 

n - 1 



n even 



odd 



(A.14) 



To calculate the T^(i)) SEK rescaled the angular sep- 
aration intervals from [t9 mm , ^max] to [—1, 1] with a variable 
transformation of the form 



2(0 - g) 
At? 



(A.15) 



with t? = (i?min+'i9max)/2, At? = tf min -if mm and calculated 
the results for t+£(x) = T+£(d). The explicit mathematical 
form of t+i(x) and i+2 (z) are shown in SEK. One can find 
the a n i coefficients in Eq. (|A.2[) by a change of variable from 
x to i3. 



The rest of the t^(x) filters are 



2n + 3 



(A.16) 



where P n +x(x) is the Legendre polynomials of order n+ 1. 
We rewrite the recursive relation for Legendre polynomials 
using Eq. (|A.15|) 



Po(x) = 1 , 



P n (x) 



1 , 




2d 


1 


At? ' 


B ' 


2n- 


1 20 


u 


At? 


n - 1 






Pn-2 


n 





Pn-l{x) 



2n- 1 1 
n ~B 



Pn-l(x) 



(A.17) 



where B = A0/(20) = (t? min - tf max )/(tf min + tax) is 
the relative interval width. Note that the first term in the 
recursive relation has the highest polynomial order, and the 



rest have lower order, respectively. We write the Legendre 
polynomials as polynomials in t? 



P n { X ) = Y,Cni # 



(A.18) 



with C n i coefficients 



Coo 
C n 

Cio 



1, 

2 

Atf ' 
1 

~B ' 
2n- 1 2 

n A^ 
n-1 

^n-2A > 



*-^n— l.i ^n— l.i— 1 



n S 



(A.19) 



where (7 n , = if i > n or i < 0. The a n i coefficients are 



2n + 3 



(A.20) 



The W„ weight functions computed the semi-analytic 
way are considerably less computationally demanding, al- 
though, formulae (|A.10[) and (|A.13|) are not stable as far 
as we investigated. The S n functions blow up for small ar- 
guments were they should reach zero and show a noisy be- 
havior. In addition, the argument for which S n starts to 
behave as it should grows with n, hence the resulting 
functions become less and less reliable for larger subscripts. 
However, that does not render them useless since calculat- 
ing the W^ ln functions from their original integral form is 
very time consuming for larger arguments were the S n func- 
tions become reliable. Since apart from S-i which needs to 
be stored once and can be loaded for further use, the time 
taken to evaluate the rest of the S n functions does not de- 
pend on their argument, which means one can in principle 
go to arbitrarily high ^-modes to calculate W^ ln (£). 

Nevertheless, in practice one does not need to go to very 
high £- modes to evaluate the integral in Eq. ([S]), and find 
the E-mode COSEBIs; since as it is evident in Fig.[IJ the 
W n functions die out rapidly at large £. By inspection of the 
plots, we deduced that the ratio of the largest peak (global 
maximum) and a peak at I « 1007rn/t? max is of around 
3-4 orders of magnitude. This property of the weight func- 
tions makes the infinite upper limit of Eq. ((5|) in practice 
manageable, i.e., the effective limits of the integral become 
finite, although they also depend on the shape of the power 
spectrum. In the present work the integrals involving W n 
functions are evaluated for a finite range of £ m in — 1 to 

^ max 

« 1007m max /t? 

max; where ^ max is the maximum num- 
ber of modes considered in the analysis of the interest. 

In the piece- wise method for calculating W n , a 
Gaussian integrat ion method (gauleg) of numerical recipes, 
iPress et al.l (|2002[) . is used for each interval considered. The 
results of these integrals are summed up as the final result. 
There is a routine in the code which finds the consecu- 
tive minima and maxima of the zeroth order Bessel func- 
tion, and puts them as the integration limits of the pieces. 
However, the low-£ values of the functions are not calculated 
in the same way, since in those regimes the oscillations of 



15 



M. Asgari, P. Schneider P. Simon: Tomographic cosmic shear analysis with COSEBIs 



T n is more important compared to Jo, and they cause nu- 
merical artifact. Instead one Gaussian integration method 
with higher accuracy parameter is used to evaluate them. 
The limit to change from one routine to the other is set by 
^thresh ~ TW^max parameter. 
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Fig. 13. Plots of figure-of-merit values with respect to the number of Log-COSEBIs modes considered for two surveys with the 
three angular ranges [l',20'] , [l',400'] , [20', 400'] . Here all the parameters except one is marginalized over. The first column 
correspond to MS (top plots) and MS+CMB (bottom plots) and the next two to LS and LS+CMB correspondingly. Eight redshift 
bins are used here. The CMB prior flattens the curves, especially in the case of MS, where WMAP7 puts tighter constrains on the 
parameters. 
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